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abstract 


tv dcvctocmem of md rasulo drained wW, a variable resolution nratcbcd-grid GCM 
for d» regional clhnana simulation mod«r= presumed. A global variable resolution stiemhed- 
grid used untie snidy bas enhanced horizomal resolution over the US. as the area of interest. 

The stretcbedgrid aptanach is an ideal tool for representing regional to global scale 
interactions. It is an alternative to the ividelv used nested grid approach introduced over a decade 
ago as a pioneering s^p in regional climate modeling. 

The major results of the study are presented for the successful stretched-grid GCM 
simulation of the anomalous climate event of the I9S8 U.S. summer drought The straightforward 
(with no updates) two month simulation is performed with 60 km regional resolution. The major 
drou„.t fields, patterns ar.d cnaracteristics such as the time averaged 500 hPa heiafats, 
precipitation, and the low level jet over the drought area, appear to be close to the verifying 
analyses ior Jie strKcned-grid simulation. Li other words, the stretched-grid GCM provides an 
ie„t uown acaj^g over the area of interest with enhanced horizontal resolution. It is also 
snown cnat u.s streunec-gnd GCM skill is sustained throughout the simulation extended to one 


vear. 


i se -evelcpec ana te^.ed in a simulation mode stretched-grid GCM is a viable tool for 
regional and subregional climate studies ar.c applications. 
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1 . Introduction 


The widely-used oesred-grid appreach represents a pioneering ^ towards 
regional climate simulations (eg. Dickinson et ai. 1989, Pieikeetal. 1992,Gionp ,990,I99J, 
Juang and Kanamhsu 1994, Caya and Laprise 1999, Pan et ai. 1999), 


Tne nested-grid approach to regional climate modeling was coming mostly fiom and was 
readily accepted and adopted by d* regional forecast or LAM (limited area model) modelere as a 
natural esrension and new application of this kind of models. The variable-resolution stretched- 
Srid (SO) approach is coming mostly from the global or GCM modeler, as a new application of 
global models in dae variable resolution frameworfc with enhanced resolution for both prognostic 
fields and boundary forcing over the area of imerest. 


1 -■* Variabie * rSO,Uti0n ^ched-grd forecast models have been first developed in the 
late /Oh. -ney are successfully used for operational short-term, 2-S-48 hour, fore-casting. A: the 
CMC (Canadian Meteorological Center), a grid point model is used since the early 90's 
(Stanifort.. -nd MitcfceJ 19/8, Scamiorth and Daley 1979, Staniforth 1995. 1997. Cote et al 
199). Ccte 1997). Independently in the late 30's the development of the variable resolution 
stretched-grid specmai model was stared at Mete^ratce donnler and Gelavn 1983, usL-tg the 
approach ceveloced by Schmid: (1977). Tne model is used (or operational short-term forecasting 
>tnc. — e ...,d-90 > (t.g. : essac and Benard 1 996). Tse first regional donate simulation was 
pertormec with the Mrjo-Frar.ce model in the mid-90's (Deque a-pd Piedeilevne 1 995). 

Cdlter variable resolution models have been developed for regional appil.-ations (e.g. 
?aeg:e .75?. Hardtxe: ..-9". McGregor and Ka czfev ! 993- c-edd-T--* -or-m r-- 

=> '° ' G ~° b) StrS:C ’ eC * gr ‘ C : - iG) * GCM is * ei "i - eve!c ?^ by me authors since the mid-90's. 

At the firs: stage of the development, the SG-zynamical core in the Heid-Suarez framework (Held 

and Suarez 1994) v, ,, . . 

::at s - V25 ce-.e:ope-c ar.c thoroughly tested in the 

~’ s * ~ ' u 1 . . v-Rac. ;tz e: a ; . . 7 9~ i. Tice dynamical 

core r.is tne ::mte-d:ne:ence dynamics identical to mat of me GEOS GCM (Suarez and Takacs 



1995,Takacs and Suarez 1996) and a Newtonian-ype physics (Held and Suarez 1994). ITje finite- 
difference dynamics has been adjusted to a variable-resolution stretched-grid and the 
computational problems related to variable resolution have been successfully resolved. In the 
follow-up study with the dynamical core, a real orographic forcing was introduced. For the case, 
the filtering technique was refined, and the down-scaling effects were confirmed. The 
experiments with the SG-dynamical core have shown that the stretched-grid approach provides 
the efficient down-scaling over the area of interest. It appears to be a viable candidate for cost- 
effective regional modeling not only for a short-term but also for medium-range and, most 

importantly, for long-term regional integrations. It was an important step towards developing a 
full diaban'c stretched-grid GCM 

Relying upon the obtained results, the ambers conduced that the global stretched-grid 
approach has to be explored as an alternative methodology to the current widely used nested grid 
approach. The following comments outline some differences of the stretched vs. nested grid 
approach. The nesred grid approach has the following advantages. The computational efficiency 
is high due to using a regional integration domain. There is the Possibility of combining the outer 
ny-rostanc GCM cr reanaiyses and an inner ncr-rydccruric model, i: is preferable to -use the 

GC. . r h> >.c to. _e i'.-er model i.e.g. Caya arc _apr.se I999,i. There aiso exists the possibility 
of using a perturbation regional model (lice the NCE? Regional Spectral Model (Juang and 
Kanamitsu 1994). 

Due to the global integration domain for variable-resolution stretched-grid models there 
-T- .A. -uUiA.a... vwi.ditiona arc no ccmputattcna. buffer around the area of interest is needed to 
con. oi me oour.ca.-y conditions problems. The global variable-resolution approach to regional 
climate modeling allows continuous/straighttbrward long-term integrations to be autonomously 

.TT. wC ] ’ t* q L it *' , * a n-a*.-* r~r~\ < 

l “- ^ a continuous or periodic updating of 

concidcni i:;ne region sroundanes. and, for ea* -. = ; ^ 

... u.xb.-w > ' — i ^ periodic re- 
initialization of conditions throughout the region :: interest. It also avoids the need to apply 



damping techniques within a computational buffer region. Both the updating and noise damping 
are required in nested-grid models to control severe computational noise arising from the 
application of lateral boundary conditions. A further advantage of variable-resolution stretched- 
gnd models is that they provide self-consistent interactions between global and regional scales of 
motion and their associated phenomena. In other words, using one variable-resolution global 
model results in a better preservation of a delicate balance between global and regional scales. 
Also, there exists the possibility of introducing a mesoscale non-hydrostatic SG-GCM. The global 
non-hydrostatic models are being developed (e.g. Semazzi et al. 1995, Qian et al. 1998, Yeh et al. 
1999), and the stretched-grid approach is a practical way to use them for regional integrations. 

The stretched-grid approach provides the possibility to use such a fine regional resolution 
that is not available otherwise for GCMs. Actually, any GCM resolution can be made at least 2-4 
imes finer for the area of interest through the SG-a PP roach. Depending on the regional resolution 
used, the computational savings are at least one order of magnitude compared to computer time 
r.eeded for the corresponding control run with a global uniform fine grid GCM. Therefore, the 
variable-resolution stretched-grid models provide a scientifically and practically attractive 
possibility of performing cost-effective regional/subregional integrations with finer resolution 
- an area o, -.rerest man are lately to be possible in the foreseeable future with fine uniform 
i.obal grid models (Cote 1997, Fox-Rabinovitz et at. 1997). 

It is noteworthy that the computational efficiency provided by the stretched-grid 
approach, is not the only reason or rational for its practical implementation. It is at least equally 
:: even more important that through this approach an efficient down-scaling is obtained 
..resencin_ iji. i.d '<ery fine regional mesoscale fields, diagnostics, and phenomena. 

All the aoove considerations explain why the variable resolution stretched grid approach 
::uid be considered as a viable alternative candidate for regional climate modeling. At the same 
V ~ *° ' S ' -■ optimal choice between using a nested or stretched-grid 



should bo candidly made for difleren. applicafioos. Also, the combinadoo of two 
approaches may appear to be an attractive optimal choice. 

Noto that according to Lindzen and Fox-Rabinovta (19S9X the coosistoucy of horizontal 
aud vertical ntsoludoo should bo pntservetL Themfme, an inmease of horizontal resolution has to 
bo accompanied by the corresponding inooaso of vortcal resolution. In dris study wo roly upon 
rather fine vortical resolution (70 layem) used for the SG-OCM experiments. 

It o worth Clarifying that the term ’down-scaling- for a dynamical model or dynamical 
down-scaling is widely used for nested-grid models driven hy a coarser resolution GCM or 
-analyses. For the nested-grid models, this term is associated with die abrupt resolution change 
between the inner model and cuter mcde, or reanalysis forcing. I, seems appropriate to adopt this 

15 a stretched grid with gradually varying global resolution. .As a 

result -he targe and medium, synoptic series are produced away them the region of interest while 
acsc.es ar„ gmdLailv introduced when approaching the area of interest (and its 

tnu_ediat. ncinity). emtciency or such a dynamical down-scaling for die SG-GCW regional 
stm-iar.on mcde is tnvasttgated in the study. 

- ht ar.cnaious .-emcr.ai climate event o: the L\S. 1933 simmer drought drew a 
-v. action Or w...u.i_oiogists and was tnoroughly investigated (e.g. Janowiak 1933. 

Rc r -.-.>sic. l 9SS. Namias 1991, Mo et ai. 1991, Oglesby and Erikson 1939, Oglesby 1991, 
tn and 3 ranstaiG r [992 A tla_s c*r 3 f t qq ^ \ -t - * * . , . 

'* ^ v ^nt was cnosen ror inter-comparison of 

nes s-*- • -rional u*oceIs. or PERCS (Project to I nrar ~onna~ a pr -■ r 

v . u 1 —v.or.pa.w K-i.onam ulimate Simulations) 

(Tamest ai. 1999 ). 


.\L-n».wei.. :r.'. estigated anomalous regional climate even: ortheL'.S. 1938 summer 

>5 rr.osen for this study for testing tire developed SG-GCM in the straightforward 

— Q ; * ■ . oce. Our ma *or - °a! is - • n - • . . 

' “ ^ i * li * fc '“ s '‘ 0i >an2Dle resolution in terms of regional 

r r — - c:r ‘S rne c;L£e sracy o: the event. 


drocont is :r 


5 [ iT. « . 
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Tie key objures of the srcdy are: developing tie SG-GCM providing an efficient 
down-scaiing over the area of inrerest wnh enhanced fine borizoptiti resolution, ami consist 

inrerections between regional and global scales; and resting the SG-GCM ability to simnlare an 
anomalous regional climate event. 

A descnpnon of the GEOS (Goddard Earth Observing System) SG-GCM is presented in 
Section 2. Section 3 and 4 are devoted to discussing the straightforward simulation results 

obtained with the SG-GCM for the 1988 summer, for global and regional fields, respectively. The 
conclusions are given in Section 5. 


2. A brief descnpnon of the stretehed-grid design and the GEOS SG-GCM 
— The stretched- grid design 

The stretehed-grid used in this study for the GEOS SG-GCM is similar to that of 

oducwv. in Siaxuforth and NGtcneli ( 1 978) and used in Fox-Rabinovitz et al.( 1 997) for the SG- 

namical core experiments (Fig. I), ine stretched grid has a uniform (latitude x longitude) fine 

r ' S0lutI0n over * e or uteres: -mch has to be a spherical rectangle. Outside the area of 

U ‘~ 5 "‘'* inter ' are 5"cressmg or stretching in both latitudinal and longitudinal 

directions as a geometric progression with the consent local stretching factor or ratio defined as 
follows: 

r j = dx/dxj, ^ 

where dx. and dx H are adjacent grid Intervals, arc j is the horizontal index. 

T:._ total globa; str etc rung tactor is defined as 

wr.ere dxnwx and dxmir are the maximal and minimal grid intervals on the sphere, respectively. 

pona=:- stretcr.ec-gnd design. the area of Interest, i.e. a uniform fine 

'••-on v val .ec wangle, can re allocated over any part of the globe such as the 

n s le L .S. _^ed in this study. Note mat if the area of interest includes the polar point 



or even the vicinity of the pole the stretched grid is to be rotated so that, for example, the polar 
point will be placed on the equator in the rotated coordinates. 

In order to keep under control undesired computational problems arised from grid 
irregularity, some important properties of the stretched-grid design have to be imposed 
(Vichnevetsky 1 987, Fox-Rabinovitz 1988, Fox-Rabinovitz et al. 1997). First, the stretching 
Should be uniform, i.e. with r ; =3D constant for all j's. Second, the stretching has to be well 
controlled or moderate in the sense that the local stretching factors should not usually deviate 
from unity by more than about 10%. This allows one to have very fine mesoscale resolution over 
the area of interest (e.g. Cote et al. 1993, 1998, Cote 1997) while allocating a significant 
percentage of the total number of global grid points within the area of interest This reduces the 
amount of computations needed over the rest of the globe. Third, to keep the overall accuracy of 
approximation ur.cer control for the moderate stretching strategy mentioned above, the maximal 

gnd haVe m ** not !ar - sr than a few (2-4) degrees to be quite close to the resolution of 

typ IC a! GCMs. This is needed for preserving the general integrity or realistic skill of global 
s.mtiat.c ..elds mat is necessary for providing consistent interactions oer.veen giobal and 
regional scales throughout 5G-GCM integrations. 

!: :J notewonhv Six: S* strictly cer.-oited SO-param«er 5 needed for the regional climate 
simdadcr. mode could he telexed for other modes of intcgmUon such as a short-term forecasting 
(Cote et _d. 1993, 1998) and data assimilation modes. 

_t.oic_ of s _ etcr.en-grid parameters depends on a particular mode; design. 
vonngura:;:n. requirements, and modes of integration. It also depends on a model's numerical 


scr.eme. 

uniform. : 
I'G-run ; 
results. T 


dU-..u,.ed ; .u '_iec m this study has the same number of grid points as the global 
— _o degree gre .. ui nut redistributed according to a SG-design (Fig. ! ). The 2x2.5 
ce.c’x ;r. section 4 as a reference run when analyzing the SG-GC.M simulation 

in:e:;St is :he 5phericai o v «r the U.S. with 60 km uniform resolution 
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and the following coordinates: from 25N to SON and 75W to 125W. For the stretched grid, the 
local stretching factors (Eq.(l)) are ~7% and -5% and the total global stretching factors (Eq.(2)) 
are -9 and -8 (that corresponds to the maximal grid intervals of approximately 5 and 6 degreesX 
for latitudes and longitudes, respectively. The stretched grid has approximately 9 times less grid 
points than that of the global uniform 60 km resolution grid. 

As an option, the spherical grid can be rotated so that the area of interest is located, for 
example, about the equator in the rotated coordinates (e.g. Takacs et al. 1994). Such a rotation is 
not necessary for the U.S. region but as it was mentioned above h makes sense for the regions 
including the pole or located in a close proximity to the pole. 


b. A brief description of the GEOS GCM 

The GEOo GCM was developed by the Data .Assimilation Office (DAO) at the NASA 
Goddard Laboratory for Atmospheres. The earliest predecessor of the GEOS GCM was 
developed in 1989 based on the ’’plug-compatible’ concepts outlined in Kalnav et al. (1989). The 
GCM was subsequently improved in 1991 (Fox-Rabinovitz. etal. 1991 , Kelrand etal. 1991). The 
Relaxed Arakaw a-Schubert cumulus convective parameterization ar.d the re-evaporation of 
railing ram are cased upon the works of Moorthi ar.d Suarez (1 992, ar.d Sue and Molod (1988). 
The .^ng-wave ana short-wave radiation is parameterized following Chou and Suarez (1994). The 
planetary boundary layer and the upper level turbulence parameterizations are based on the level 
2.5 closure model of Kelrand ar.d Labraga ( 1988 ) and Kelrand er ai. ( 1991 ). The orographic 
<_.ag r a.a...e:enza:;on follows Znou et a!. ( 1995 ). The model physics is updated 
.-.-..t i..w_,..c.es lai.gn.g from ever, two cynamics time steps for turbulence and gravity 
wave crag, three dynamics time steps for moist processes (convection and large-scale 

- . ). t„ i.cl, .or s,.ort-wave ar.d three r.curs for long-wave radiation. All model 

pnystes updates are prorated ar.d applied at every time step. 


9 



The momentum equations used in the GEOS GCM are written in the "vector invariant" 
.orm, as in Sadoumey (19/5) and Arakawa and Lamb (1981), to facilitate the derivation of the 
energy and potential enstrophy conserving 4th-order differencing scheme. The thermodynamic 
(potential temperature) and moisture (specific humidity) equations are written in flux form to 
meditate potential temperature and moisture conservatiomThe Arakawa C-grid is used for 
Horizontal approximation. For vertical approximation an unstaggered Lorenz (I960) grid in 
generalized sigma coordinates, is used. The vertical differencing scheme is described in Arakawa 
and Suarez (1983). The time integration is done with the economical explicit scheme (Ames 
1969, Schuman 1971, Brown and Campana 1978, Fox-Rabinovitz 1974), based on a leap-frog 
scheme with time averaged pressi^e gradient The scheme allows one to use approximately twice 


T'"‘" dm.- sieps than chose of the leap-fros scheme. 


complete description ot the fourth order 


finite-difference scheme used in the dvnamical 


core, can be round in Suarez and Takacs (1995). 


Since the early 90s. foe GuOS GCM is routinely run with 2 x 2.5 degree horizontal 


resolution and 70 layers in foe vertical 


covering foe entire tropospheric and stratospheric domain 


~er.vsen the sl rface «* the 0.0! h?a leva!. The GEOS GCM ■•vas developed and esed for 


:::cucing global rear.alvses ford: 
23 Loose or the A\fT? - Acmes 


ate applications. It is also used for long-term simulations 
'.eric Model [nter-comparison Project). 


c. Configurations of foe GEOS SC--GCM 

A: me current stage of foe SG-C-CM development, a stretched grid is fully implemented 

.a....c_, . -ai c ro graph;, as described by Fox-Rabinovitz et al. (1997). 
-..mg wim introducing a wretched-grid coordinate array, some interpolation techniques 
depending on a grid definition were modified accordingly. The model filtering procedure has also 
rn ° U - :3C : ° 2CCC -" : rcr ** changes required for a stretched grid by Fox-Rabinovitz et al. 
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Tie prognostic vatiab.es (wind consponeoa, reospentinre, motauns anti turfite* preK me, 
tor the dynamical mode, state am updated and stored on the model dynamics stretched grid. The 
diabetic tendencies are updated ax the appropriate physicalcompumtionai time-sca.es on their 
own physical (imennediare uniform) grid. Tien they are imerpolaed, pmraxed per time smp and 

applied a. eve^r time smp » the mode, dynamics stietched grid. Therefore, me whoie inm^tion 
history resides effectively on the stretched grid. 

Such an approach is justified by tire assumption that model physics and dynamics can 
have ddfermt temporal and spatial resolution. This subject has been discussed by Lander and 
Hoskins 0997). For a special model, drey advocate using coarser resolution for model physics 

titan for model dynamics, and further conclude that similar considerations also apply m finite- 
v^rrersnce and firute-eierrienc modeis. 

A. this tim >*age of the sC-GCM development Lmplementing model physics on an 
intermediate uniform resolution grid allows us to avoid some potential complications that may 
-.se from calculating all model physics parameterization* on a stretched grid. It was verified that 
-ombmanon ol the Su-rr.ede! dynamics and intermediate uniform grid model physics, 

- sico ( ices . me rir.er scale patterns produced by the model dynamics on 

- attw.v.ned gnd tig.2 snows the instantaneous regional patterns of specific humidity at 850 hPa 
obtained on the stretched-grid (the upper panel) and interpolated onto the uniform intermediate 
caution grid used for model physics calculations!^ bottom panel).Both patterns show similar 
gradients arc other features. This similarity takes place because the model dynamics and 
orographic forcing arc the whole integration histoty resides effectively on the stretched gnd. 

At me next stage of the model development, the model physics parameterization* will be 
gradually (cre-by-or.ej implemented on a stretched grid. 

Ore granny is calculated directly on a stretched-grid by averaging, within a gnd box. the 
* - ... 0 p artace e.e-.at:on dataset available from the National Center for 

ye-... R. .N'CARj. i r.e gnd box averaged orography is passed through a 
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Unczos(I966) fiber m bo* dimensions which removes the smallest scales while inhibiting Gibbs 
phenomena (Takacs et aL 1994). 

The ««d»d grid simulanon is compared below to that of a reference run performed 
die global uniform 2 x 2.0 degree resolution grid with the same amount of grid points as that 
of a snatched grid. A strict ultimate validation apptoach used in this study is comparing the 
simulated fields against data assimilation, or reanalysis products. It allows us to calculate the 
mean errors or biases as the deviations of simulated fields horn the conesponding verifying 


analyses. 


'■ G! ° bal SeaSOnaI fidds for 1116 summer simulation with the SG-GCM 

we star: analyzing ±e results of tine SG-C-CM simulation with presenting the global model 
Saids for the 1933 summer season (fee- July-August, or JrA). In this sectron, we will focus on 
discussing the global Integrity of simulated fields although the SG-GCM has been developed with 
the emphasis on regional climate sundadon. In the context of regional ciimate modeling based on 

—c.. _.e ..wgrfty of global fields :s ar. enporaht property of an SG-GCM 

*"" U " ' ** " oul !n — previous Section -using a variarie-resolution SG-GCM 

■* : ° ?r ° V:de -P ro - Dnate --spresentaticn ofgiccai large scales and their consistent 

interactions with regional tnesoscales. 

The SG-GCM integration started at 00Z May 1 5, 1 933 and confirmed through August 3 1 . 

: 9i3. The mode! progr.csric and diagnostic fields a-, stored every 6 noun and used for producing 
e-a'. eragec ::e ; as. The SG-GCM simulation is run with real surface boundary conditions 

m °nthlv ...ean analyses or sou moisture and weekly mean analyses of 
~ 1 ' 5ea SUrfaW tem - ce — >• Snow * ~- d sea ar; ««*. The GEOS DAS reanalysis products 

0tTn “ C : degr ” ’-soutior.. are used for alidation of the SG- 
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One of the major goals of the SG-GCM simulation presented in this and the next section 
15 10 WhCtfaer ** results obtained earlier with the SG-dynamical core in the 

Held-Suarez framewcrk (Fox-Rabinovitz et al. 1 997), will be also obtained for the full drabatic 
SG-GCM. Namely, two major results have to be verified for the SG-GCM simulation: the 

efficiency of the down-scaling over the area of interest, and the overall 'integrity" of global 
simulated fields. 

Let us start discussing the results of the SG-GCM simulation from comparing the 
obtained seasonal zonal mean vertical distributions for dependent variables against those of the 
verifying GEOS DAS reanalyses produced for the 1988 summer. The I1A wind component 
distributions ate show, in Fig.3. and those of temperature and telarive humidify in Fig.4. 

The SG-GCM simulation distribuncns appear to be quite similar to those of the verifying 
reanalyses. The positron and strength of the jets in both hemisphems are quite well reproduced in 
u.e SG GCM stmulatrcn (rtg.ja and b).The Southern Hemisphere jet core (the 40 m/s isoline) is 
stronger oy 10 m/s, laner, and shifted equamraard in the SG-GCM simulation (Figsja and 5a) 
compared to shat of or -he reanalyses (rig.4c). The positive bias (or the deviation of the SG-GCM 
S"..uiat...n corn the veaytsj GEOS analyses) of -0 ms in the Southern Hemisphere tropics 
-outu _s aub— 00 c-a taver (Tig. 5a). is associated with the equatorw a-d shift of he jet. The 
larger bias takes place above 100 hPa. All he differences result from coamer resolution within the 
a..tched grid .or the s-buthem Hemisphere. L, the better resolved within the stretched grid 


Northern Hemisphere, the je: position and strength <Tig.3a) are dose to those of the reanalyses 

^ ^malle*. -r.ar. : m-s. except tor the jet core where it is about 


‘ “ e rn ^‘ wIOna ’ v!nci distri b u tions :cr both the SG-GCM simulation and the GEOS 


~ -*ose to eacn cuter f Figs. 3c and d). Th 


e similarity' taxes place in the 


2 - ::er *5- - - - s and in both the latitudinal 


; posiuon in tr.e tropics anc in the vertical 


^re also close for both the 
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and mamtiyai a Scida. The bias in dm Norton, Hemisphere is small. less than 0.5 m/s. ,. 

' S * ^ blSSer ' ab0UI *' °' 5 mS - ” *' tr °P ics for the 300-100 hPa layer, especially in the 

Somhon Hemisphere, and widbm the somhem poiar cap for the lower tiopospheric ,evels (Fig 
5b). 

The cold tropical tempemrure core as dte 100 hPa levei. the cold polar tempemtin* core 
above 200 hPa in the Southern Hemisphem, die posidon of the maximum smface tempemmre m 
the Northern Hemisphere tiopics (the 295 K de^ee isoline, as wel, as die typioai gradients, show 
lose partem stmilariiy for the temperature distributions of the SG-GCM simulation 
and the verifying reaves (Fig4a and b). The temper bias (Fig.Jc) is less than 3(K for die 
ennre verdcai-larmrcinal domain, with the exception of dm positive (warm) bias of 6(K around 
-e 100 h?a level m .... tropics. The larger warm bias up to 12(K is coniined to the Southern 
Hemisphere coiar cap and located above the .00 hPa level. Same hind of the overall pattern 
similarity is obtained :or relative bumidiry including the positions of the maxima near the surface 
especially Ln the Norhem Hemisphere (Fig.4c and d). Toe relative humidity bias (Fig.5d) is quite 
small, mostly w~fc '.0-10%, for the entire Norhem Hemisphere except tor the poiar cap around 

"* ‘"" 2 ■ SV * J ' r ‘ e fci “ :s * to lar *" v to 30’-,. around 50(S at 200 hPa and up to - 

J/ * near : "' e Sc --~ around 500 hPa. 

Tr.e norzomal distributions of relative humidity at 850 hPa for the SI A SG-GCM 
simulation, die '--nrying reanaiysis and their difference or bias, are presented in Fig.6. The 

reanaiysis fields (Fig. 6a and b) show the close similarity with each other, 

I> ... ...e Nomem rtermsphere. Tne bias (Fig.fic) over the area of interest and its vicinity 
.s small, mostly : 5* and below indicating that the efficient down-scaling paxes pl a ce for ^ ^ 

• •~- 0 ^ u.w pa_ems over me southern Hemisphere are also quite similar (Fig.6a and b) the bias 

! ~ : X. 6c i is ' --'-a* ^ ^ \* • r w . 

**“■** Non^em Hemisphere. 

l5Sw " i 3r ‘°' V tna: the overa11 '-xgny of global fields is preserved. 

>01 red w:-nn the stretched grid Northern Hemisphere. The larger bias or 
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dK devi*™ of d. SG-GCM timulntion flelds fora ^ oftfe fc ^ 

Scothom Hanisphert does not seetn to direct* negative* affect the resuta over the are. of 

also have to realize that some deviations are not relaxed to using the SG-GCM but are 
inherent to the basic uniform grid GEOS GCM_ 


4. Regional fields for the 1988 U.S. summer drought simulation with the SG-GCM 

in this section, we wfl] discuss the regional fields representing the drought. For analyzing the 
down-scaling effictency over the area of interest we compare the SG-GCM and the 2 x 15 degree 
reference UG integrations against the GEOS DAS verifying reanalyses products. For the 
validation purposes, we wifi mostly use the June 1 988 monthly mean characteristics and also 
these or the two month drought period from May 15 to July 15, 1988. 

t L.e extent of _e anomalous drought is shown by the anomaly distributions for June 1988 
for me 500 hPa height and total precipitation fields (Fig.7). The anomalies are calculated as the 
differences oerween the monthly mean obtained from the GEOS DAS reanalyses and the June 
climatology. The 196X995 climatology is used for calculating the precipitation anomaly 
'Tlg.-a, arm that ofl9*X?95 for the 500 h?a height anomaly (Fig.7b). The NCE? gauge 
pre_.piui_.cn ;or June Ir?3 show in Fig. 7c. 

The sow* positive er.oma* for the 500 hP. heights over the U.S. and Canada reaching 
is accompam*,. dv two negative anomalies, one on the east coast and another just off 
-a. -v. voti. r.a * cp to oQ m (Fig. 7b). The strong negative precipitation anomaly cr 
~ ~ ' a ^ ‘° C 0Vir uUie Mississippi River basin reaching up to -2.5 mm/dav 

(F i= . _). eganve arommy over the Mississippi River is a strong anomalous regional climate 
-w 5 - Xouelewic .*33. Janowiak 1983). The more moderate positive rainfall anomaly 

wa - ?!acs over :he south-central par: of the C.S. The Midwestern 

crc-gr.t area .vith orecm.mr less i ■ c . _ 

-vis . . i..m t-i_. .5 jno’-'r. in Fiz. /c. 
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a. The 500 hPa height fields 

Let us compare the 500 hPa height regional pattern produced by the SG-GCM simulation 
against the GEOS DAS verifying analy ses. The two months mean field obtained from the GEOS 
DAS (FIg.8d) shows that the anomalous drought event is characterized by a strong ridge over the 
central part of the U.S. and Canada and by two strong troughs located on the east and west coasts. 
Sucn a pattern is consistent with the strong anomaly shown in Fig. 7b. 

The two months mean oOO hPa height field obtained from the SG-GCM simulation is 
presented in Fig.8a. It exhibits the major partem similarity to the GEOS DAS results (Fig.Sd). 
Namely, the SG-GCM simulation (Fig.Sa) contains the appropriately positioned ridge and two 
troughs although the features are less pronounced than those of the verifying analyses (Fig.Sd). It 
seems appropriate :o emphasize da: this is a positive result for the straightforward two months 
sim -teucn. To assess it quanctarively, the bias or deviation of the SG-GCM simuiadon (Fig.8a) 
from me veruying analyses (Fig.Sd) is shown in Fig.Sb. The maximum negative bias values 
reacting -60 m to -30 m are confined to the relatively small area immediately adjacent to the 

. ..orthem „oui.da._. of me area or interest. Over the remaining Iarze part of the 
ar__ -. zte oias va:ues are much smaller and are of the different signs. The bias here is 

most.;, within me -20 m to 20 m range. The bias distribution is consistent with but definitely 
smaller than the June 1988 anomaly (Fig. 7b). 

Tne etas .or the reference run or the UG-simulation is presented in Fig.8c. It is negative 
over me entire area of interest and is reaching -i 10 m. For the large area south of this minimum, 

m - * ,J - •• * lomg d ov. to -4v m to -20 m over the 

....... area Oi i es^ me ov era.. LG-simuiauon bias (Fig. 8c; is approximately 

r ’ V ° ‘~- er rnagnituGe :han Lha: of -“- e SG-GCM simulation (Fig.Sb j and is comparable 

vi! " ” : '- r ‘ ; « 3 ! " omal > All :r,; *, SG-GCM sim.Atior darimtdy 

d0S!.- :: As -trir- :.-.j analysis f* UO^fewtoSw for 500 h?a 
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W= compansd the 500 bPa regional voting analyses ^ ^ ^ ^ 

^ ^ GE0S S °- DAS — - - - ■*- «**a They „ „ 
*e * regiona, wi* * same ,^ 0 „ of ^ ^ ^ ^ ^ ^ 

" FiS - 84 ^ bia » P^ «*« ^ finer msotadon analyses am of .he same 
magnitude as shown in Fi» 8b. 


b. Precipitation 

°ta«n S * «** F-W-ta P*te „ evidendy d* most taportant ^ of . fc 
dro^n stmulanon. I. is dsc da. most diffiod, fieid „ simidae due „ die convexity of*, 
process. Therefore. when analysing predpimaon dishibimons. we have ,o concernm* on the 

^ ^ 68 ^ •«— * «— ** straightforward w 

^ 1988 — *«*»*. *>* - biases fo. done 198S „ 

i-r.semea to F.g.9. The precipiahon partem produced by die GEOS DAS (Fig.9d) shows dearly 

dh*v*- area for the Mississippi River basin. The minimum (less than I mm/d i •• ■ 

u man 1 nzm/day) precipitation 

*- ceased in thenonr.em part of the basin The -d-c— ? nr»r' ■ • 

. ^ 1 au "- ? r -cip!tat:on occurs over the major 

-ew me lh tains. .A stronger creci-irar ; cn ^ ^ ^ . 

~ ,c ^ a >- ; s octamed over the south- 

1-17 K ^TT ^ C - 5 - ^ PreCi?itati0n *«“"* » consistent 

‘y i -a.on ancta aly ana preapitation panem shown above (Fig. 7a and 7c). It is 
nctev. onhv llthoi'an n ^ c - . 

' ^ DAS preCIp,laCl0n <***> *"* -me deviations/hiases from *e 

NCc? gauge precipitation (Fig.7c) thev are loc-^ n - , - . 

- *<*•— *™sd, OLuiide or me drought area. They do 

“C" seem to afreet si^-im •- 

' uLl -iscussion or the nG-C-CV dm.— -r -• • • 

urou^nt simulation. 

‘ ^ "° n '* a " W '‘ °°~ ir ' ea rr ° m ** SG * GCM simulation (Fig.9a) shows the 
"“ UC "' preC ’ F!tar:0n - beiow 3 m ^ay. over the major drought area, the nonhem 

5 ^^3S:5S:nnj 

‘ Wie rna ior part or me Rocky mountains, precipitation is 

* ? ” 4 ° ^ ~ keS ?;ace over *e south-central, eastern and 
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pmufd* as. AH mate *e preciptauou patjern obtained from fc SG . 
GCM simulation, quasi-reaiistic. 


Tr.e : 


P^oprancn biases or deviations calculated against the verifying diagnostics 
P™hiced by the GEOS DaS, for the SG-GCM simuteoo and dm referance UG-simufction, are 
»9b and Rg.9c, correspondingly. The SG-GCM simulation bias (Fig.9b) is within 
*e 2-3 mmfday tange over the draught area and is definftely smailer than that of the tefetence 
UG-simuiarion tanging there from 3 tnm/dav to 6 mm/day (Fig.9c). Such a significant bias 
station for the SG-GCM simulation precipitation compared to that ofthe reference UG- 

simulation. is an important indication of an efficient downscaling obtained for tins anomalous 
drought event through the SG-approach. 


■ - LJ) makes a significant contribution to the drought formation and 
presentation is important for a successful regional climate simulation 


c. Low level jet 

The low level : e 
sustenance. i-5 axi e^u a ie 
or the anomalous ever.i 

" " '" iar ' ''' inC V " :crs :or June i9ss representing the LLJ for the SG- 

C-CN. -cne me rererer.ee L G-simuIaticn. and for the verifying GEOS DAS reanalyses, are 
sr.cwn tn F ;g. ! 0. According to the verifying analyses (Fig. 10c), the southeastern flow 
approaching Texas mom rite Gulf of Mexico gradually turns clockwise. It becomes a 
predominantly southern few east ofthe Rocky mountains or over the western pan ofthe 
•r..‘ -y- • - - .... n. er me central and eastern pan of the basic, the flow is weak. 

- well reprccuced in the SG-GCM simulation (Fig. 10a). The 
somheastem. flow over tire Gulf of Mexico turns clockwise and becomes the southern flow over 

T ““ "* ** ** ~ “irccticr. and intensity similar ,o those of 

:-.er the central and eastern pan of the 

" DaiI “ ' 15 ^5imi. a: to that ofthe verifying analyses Fig. 10c). Even 
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outside the area of interest for the immediate vicinity of the U.S„ the SG-GCM simulation wind 

velocity is quite similar to that of the verifying analyses, with the exception of the flow over the 
central and eastern Canada 


However, (he LU panem for *e reference UG-simulation (Fig. 10b) is completely 
difierem Sec that of the verifying analyses (Fig 1 0c) ovor *u draught area. Tbete is no southern 
Bow located east of the Rocky mountains. Instead, die southeastern flow coming from the Gulf of 
Mexico to Texas utras westward. As the resole the strong westerlies Bow throughout the entira 
draught area and further east. Such a LU panem for the teference UG-simulation (Fig 10b) 
appears to be not as realistic as that of the SG-GCM simulation (Fig.IOa). 


Let us consider cow the diurnal LU variations. The nocturnal LU is characterized by the 

maximum "-ds taking place around 3-5 a.m. • local time). To show the diurnal variations, the 

monthly mean wind veemrs are calculated for 10Z and 12Z. These OOZ and 12Z monthly means 

are overlapped tor comparison purposes, and presented in Fig.l I. The overlapped wind vectors 

are ’V-shaped for the southeastern flow comings from the Gulf of Mexico to the western pan of 

me drougm area. The tight ^"-branch corresponds to the I2Z monthly means and the left V- 

— c.t .0 ,Z aiu.... ...eans. For the area :: me drought I2Z corresponds to -5 a.m. local 

mme for which the noerrma: LU winds are supposed to be smonger than those ofCOZ or ~ 5 p.m. 
local rime. 

The wind vectors obtained from the verifying analyses (Fig. 1 lb) for the LLJ east of the 
Rocky mountains, are rotated clockwise from the left OOZ "v ’-branch to me tight :2Z V’-branch. 

- r.e r.gn.t 1 2Z '’- ’-branches are mostly longer man. the OOZ ones that is appropriate tor the 
.xv.uma: — . ; ms major rearms as .veil as the above clockwise rotation are well reproduced in 
. airmJatio.. . ig. . i a). The whole panem for the oG-GCM simulation is quite similar 
:o mat o: me vetitying analyses (Fig. Ilaandb. The pattern and the above features over the area 


reference LG-simulatic: 


- 1 - ik.OLvn;. -v::n :r.e cc...:na::ng -ves:er!ies over 
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the drought ana (Fig. 10b), are significantly different from those of the verifying analyses 
(Fig- lib). 

One of 4e plausible explanations for the differences between the LU patterns obtained 
for the SG-GCM simulation and the UG-simulation, could be the feet that orography is better 
resolved for the stretched grid than for the UG within the area of interest As the result, the SG- 
mountains are producing stronger blocking effects for the flow coming from the Rocky 
mountains into the drought area. The weaker blocking effects for the UG-simulation result in 
strong westerlies developing over the drought area that is not supported by the verifying analyses. 

Similar results consistently favorable to the SG-GCM simulation vs. the UG-simulation 
are obtained for other model fields and diagnostics. 


d. extended simuiaticn and dovve-scaiLns effects 

i ne bG-GCM simulation was extended to one year to verify whether the positive down- 
scaling effects discussed above will persist in longer integrations. This is an important issue for 

potential applications trie 5C-approach to long-term climate and climate change regional 
assessments. 

i i.e exter.de c simulation results are presented for the winter season or DJF (December 
I9S3, January and Fefcraiy 1 939;. Note that the qualitatively similar results are obtained for the 
fall of 1938 and the spring of 1989. The overall integrity of global fields for DJF is similar to that 

P l '’ Se ‘""* ^ ^ e -ill focus here on discussing the regional down-scaling aspects 

ror DJF. 

D.. -L. e SG-GCM simulation and the reference UG-run are compared 

against me ver tying GEOS DAS analyses in Fig. 12 . The LLJ for the SG-GCM simulation 
<r:g. 12a; anc me verm r.g analyses i. Fig. 1 2c > are close to each other. Both patterns contain 
i.. ...g j.ju.i.em :t mmg mom Gulf of Mexico. Tne UC-run has failed to rer resent this 
.....or.a... feature and sr.:w S only me westerlies prevailing over the region east of the Rockv 


20 



Mountains (Fig. 12b). Tin prognostic variable biases for difierent model levels are consistently 
smaller for the major part of tbe region by approximately 30%-10% for *e SG-GCM simulation 
comparod to those of dm UC-ron. Tbe biases for 350 hPa tempcmtums are presented in Fig. 13. 

Tbe SG-GCM simulation has smaller biases (Fig., 3a) than those of tbe UG-run (Fig. ,3b) for dm 
region east of the Rocky Mountains. 

To continue tbe bias discussion we ptesent the global zonal mean DJF biases calculated 
against verifying analyses, for tbe 925 hPa meridional winds and 350 hPa specific humidity for 
the SG-GCM simulation and the UG-run (Fig.I4). Over the tegion of interest (125W to 75W) the 
SG-GCM simulation biases am predominantly smaller than those of the UG-run. Outside of the 
area of .merest, the UG-ron biases am mostly smaller. The SG-biases in the S. Hemisphere 
depend on coarser resolution within the stretched grid compared to those of the bear resolved N. 

.he.^ I L s -o - e -.e significant dependence or. local, regional resolution for the SG- 
C-CM simuianoc. 

The legitimate puesuon to ask concerning me betier skill SG-GCM simulation whether 

me efficient dowit-scaiing to r.esoscales will be seen in seasonal mean Sei is. In other words, will 
-i.Tie averanr ~ 3 _ _ _ 

" = structures prccucec La Lae course o: integntien or they 

" ?er5,StSn: ir ‘ CUgr ‘ tC ~~ 5ea50nd “ T - e 5£a!iS? 525 h?a DJF veracity fields for the SG- 

imUiar:Cr ‘ * u ‘ € LG * rjn 212 pr2Sented u " F ^- 15 - Tne mesoscaie structures are produced 
for the SG-GCM simojation over die Appalachians and die ms. of die region east of die Rocky 
...ountains (F.g. I ;ai. Omy synoptic scale features are produced for the CG-ron .'Fie. ' 5b). 

.< s.<... .ur u._ ou-GCM simulation e ccr.ciude Laat the 

" “‘ a " :n:rc ^- c: "= sriianced regional resolution resulted in 

cer.er regional simulation resuits. 
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5. Conclusions 


TTie developed SG-GCM is the new tool for producing enhanced resolution regional (not 
global) climate simulations. The hypothesis verified in this paper is whether the potentially 
posmve impact ot tmemmhanced resolution provided for prognostic fields within a stretched grid 
over the region or interest, can counteract the potentially negative impact of a relatively poorly 
resol ed flow coming mo the region. Actually, a 'working compromise* not a 'perfect solution" 
is sought for the region of interest when applying the SG-approach to regional climate modeling. 
The results obtained in the study show that the developed SG-GCM provides such a 'working 
compromise so that the SG-GCM regional fields are closer to the verifying analyses than those 
of the reference UG-run. Also, the efficient down-scaling to mesoscales is obtained over the 
region of interest for the SG-GCM simulation. 


The following conclusions are arrived at as the result of the study. 

I . The variable resolution stretched-grid modeling approach is introduced to the GEOS 
GCM with ^difference dynamics for the regional climate simulation mode. The 
-iraore resoumcr. sremhed-grid version of the fully diabatic GCM is developed and 
successfully tered in the simulation mode. 

2- Tice straightforward continuous (with no updates) stretched-grid GCM simulation with 60 
!<m regional resolution is performed for the 1983 summer. The seasonal global fields 
show the overall global integrity that is necessary for an adequate representation of 
vons.i.en: m.te.-acnons between global and regional scales. 

i>.~ tu..wC-i..u GCM summation or the anomalous climate event of the 1988 U.S. 
SUITU„». ioUi... shows me erticient down-scaling over the region of interest. The SG- 

7531:1:3 “ e c!oser t0 --ose of the verifying analyses than the 

resul-u me reference experiment with the global 1 x 2.5 degree uni term grid with the 
same -r.cur.: c: gr.o points as for the stretched sf.z. 
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4. the major drought characteristics such as precipitation, the 500 hPa height pattern, and 
the LLJ, are successfully represented in the stretehed-grid GCM simulation, 
o. Similar positive down-scaling effects a ire obtained for the emended to one year 

simulation. It shows that the SG-approach is valid for loog-tenn regional climme 
simulation and related applications. 

6. The obtained results show that the snetched-grid approach is a viable candidate for 
regional climate simulation studies and applications. 
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Figure captions 

F ‘g- 1 The stretched grid with the area of interest over the U.S. 

Fig.2 Eelds of 850 ^ specific ^ (spHu _ 85o) w ^ 

SmU “°° Whh ' “ ' ta ra0to ' M * — dynamics and ft) the same field 

"“ P ° tetd ^ “ i “ OT, ' diatt “ nif0 ™ «*— grid used for Ccfoaring model physics. The 
contour interval is 1 g^fcg. 

F.g.3 Scasonai (JJA) wind component zonal-vernca, distributions for (a) U-wtod for dtc SG- 

OCM simulation; ft, U-wind for the OEOS DAS reuses; ft, V-wtad for ft. SG-sfotuiarion; 

(d, V-wind for ft. OEOS DAS ^analyses. The comotm m^vais are 5 m/s for (a, and ft* and 0.5 
m/s for (c) and (d). 

F'g.4 Same as in Fig.3 bn, for temperature, (a) and (b), and relarive humidity. ( c , and (d). The 
contour intervals are 5 K degrees for (a) and ft) and 10% for (c) and (d). 

Hg.5 Seasonal (JJA, aonai-vetrica, bias disrtbunous or foe deviations of the SG-OCM 
sunulanon nelds 6om foe foe GEOS DAS reanaiyses (both shown in Figs. 3 and 4) for: (a) U- 

•vmd, ft) V -Wind; (c) temperature; and (d) relative humidity. The contour intervals are 5 m/s for 

(a) , 0..' - s ror (b): 3(K for and 5% for id). 

F '°' 6 ' 5i0 " al (JJA) norzontai distributions of relative humidin.- at 850 h?a for: (a, foe SG- 

GCb, simulation; ft, foe GEOS DAS reuses; and ft, foe bias or differences between (a, and 

(b) . The contour intervals are 10%. 

' 2> ?r2CiPltati0n r ° : ' 500 hei ^ Amalies, and (c) NCHP gauge precipitation for June 
-entour in... .al> aeO.o mtn/day ror (a). 10 m tor ft), md I mm/day for(c). 

F ‘ g 8 *** M# m0 " ,h (May 1 5 •» Iuly U, 1983) 500 hPa heights for: <„ foe SG-GCM 

vunulaucp: to) the SG-GCM simulation bias: (c) foe reference UG-simulation bias: (d) foe 

Dm. C--OS DAS Italy ses Tne contour interval is -to m for (a, and (d): and 20 m for (b) and 
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Fig.9 Same as m Fig.8 bet for monthly mean (June 1988) precipitation. The contour interval is 
I mm/day. 

Fig. 1 0 The monthly mean (June 1 988) wind vectors for me 925 hPa level or LU for (a) the SG- 
GCM simulation; (b) the reference UG-simuiation; (c) the verifying GEOS DAS reanaiyses. 

Rg. 1 1 The 00Z and 12Z monthly mean (June 1988) wind vectors (overlapped) for the 925 hPa 
Wei for (a) the SG-GCM simulation; (b) the verifying GEOS DAS reanalyses. 

Fig-12. The seasonal DJF (December 1988-January-Februaxy 1989) mean wind vectors for the 

92o hPa level or LU for (a) the extended SG-GCM simulation; (b) the extended UG-simulation; 
(c) the verifying GEOS DAS reanalyses. 

Fig. 1 3. The seasonal DJF 850 hPa Maparmare biasas calculated agatat the GEOS DAS 
vaAA.Ing analyses :br (a) the axan&d SG-GCM simulation; and fb) the extended UG- 


ydnulaticn. The contour interval is 0.5 decree K 

Fig. 1 4. me DJF zonal mean biases calculated against the GEOS DAS verifying analyses for the 
extended SG-GCM simulation (open circles) and the extended LG-run (dark circles) for: (a) the 

9 "' hPa rnendl ° r * a: WKd ^ ^ *-- e 850 specific humidity (SPHL'-850). Tr.e units 


" SeC 2rC fcr -~ e clonal wind and specific humidity. correscondmzly. 

* i 5e ~° r ‘ al ° jr 9:5 h?a VCKicit >' for: & ** extended SG-GCM simulation; and b) the 
extended LG-simuiaticn. The contour interval is 10' 5 s' 1 . 
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b) 500 mb Geop. Anomaly - June 1988 
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